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ABSTRACT 

The role of acoustics in flame/vortex interactions is examined via asymptotic analysis and numeri- 
cal simulation. The model consists of a one-step, irreversible Arrhenius reaction between initially 
unmixed species occupying adjacent half-planes which are allowed to mix and react by convec- 
tion and diffusion in the presence of an acoustic field or a time- varying pressure field of small 
amplitude. The main emphasis is on the influence of the acoustics on the ignition time and flame 
structure as a function of vortex Reynolds number and initial temperature differences of the reac- 
tants. 
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1. INTRODUCTION 


Acoustic waves play a pervasive role in compressible flows even in the low Mach number limit. They 
are a necessary concomitant of nonlinear interactions such as those present in turbulent flows. For instance, 
the interaction of any disturbance, whether it is of the entropy or the vorticity type, with steep gradients as 
in shocked flows [e.g., McKenzie and Westphal(1968), Zang, Hussaini and Bushnell (1984)] generates 
acoustic waves. Chemical reactions enhance them [Toong, et al. (1974), Roberts (1978), Clarke (1985), and 
Jackson, Hussaini and Ribner (1993)]. Conversely, acoustic waves can engender instability in laminar shear 
flows which may eventually suffer transition to turbulence [Goldstein and Hultgren (1989)]. In reactive 
flows, they play a similar critical role. It is now well established that they can significantly affect the stabil- 
ity of flames [McIntosh (1986)], the transition to detonation [Urtiew and Oppenheim (1966), Lee (1977), 
Shepherd and Lee (1992)], etc. Oran and Gardner (1985) provide a fine review of combustion-acoustics 
interactions with an emphasis on the physics rather than the mathematical analysis and techniques. 

In the present study, the role of acoustics in the flame/vortex interaction process is investigated. In 
order to provide a proper perspective, the literature on the flame/vortex interaction problem is briefly 
reviewed here. This problem was first formulated by Marble (1985) to model one of the basic physical 
mechanisms underlying the complex processes of turbulent diffusion flames or combustion in vortex- 
dominated flows. The model consists of a diffusion flame with fast chemical kinetics (along the horizontal 
axis separating two reactants occupying the upper and lower half-planes) which is distorted by a vortex with 
its center at the origin. The theoretical results have established that the flame sheet is rolled up into a spiral 
around the vortex on the convection time scale forming a reacted core, and then spreads across the spirals 
on the diffusion time scale; the growth rate of the reacted viscous core obeys a similarity law as well as the 
reactant-consumption rate which is independent of time. These physical features are confirmed by numeri- 
cal simulations [Laverdant and Candel (1988), Rehm, Baum, Lozier and Aronson (1989), and Norton 
(1983)] implying, thereby, that the inclusion of shear (which was neglected by Marble) or finite-rate chemis- 
try does not alter the qualitative picture. The premixed combustion in a vortex also displays similar features 
[Peters and Williams (1988)]. 

The aforementioned studies assume the simultaneous existence of a flame and a vortex and the subse- 
quent evolution. Thus, they preclude the ignition regime which is the main thrust of Macaraeg, Jackson and 
Hussaini’s work (1992). Their work is based on the assumption of a constant density, one-step Arrhenius 
reaction between the fuel and the oxidizer, occupying adjacent half spaces, in the presence of a vortex with 
its center on the axis of demarcation between the reactants. Their study focuses on the ignition time, loca- 
tion and the flame structure as a function of the vortex Reynolds number and the initial temperature 
differences of the reactants. In the absence of the vortex, the problem, of course, reduces to the classical 
case of Linan and Crespo (1976) who analyzed the continuous temporal evolution of such a configuration 
from nearly frozen flow to near equilibrium flow. The reactants mix by diffusion until, at some finite time, a 
thermal explosion occurs at a well defined location, and ignition takes place. After ignition, a pair of 
deflagration waves (or premixed flamelets) emerge according to classical thermal explosion theory. These 
premixed flamelets are quite weak in that the temperature rise associated with them is small, and they exist 
only until all of the deficient reactant is consumed. Just beyond the deflagration waves, a diffusion flame 
regime exists where the mixing process is governed by diffusion in the direction normal to the flame. As 
time increases, the diffusion flame approaches a flame sheet We note that the existence of a well defined 
ignition point and the premixed flamelets depends critically on the relative magnitudes of the two initial 
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temperatures to that of the adiabatic diffusion flame temperature. That is, if the adiabatic diffusion flame 
temperature is greater than either initial temperature, a well defined ignition point always occurs, followed 
by the premixed flamelets. On the other hand, if the adiabatic diffusion flame temperature is between the 
two initial temperatures, there is no well defined ignition point, and a single premixed flame merges 
smoothly into the diffusion flame. This configuration of a triple-flame (or tribrachial flame) appears to be of 
such a fundamental nature that it also exists in steady supersonic reacting laminar mixing layers [e.g., Jack- 
son and Hussaini (1988), Grosch and Jackson (1991)] and flames propagating into a nonuniform mixture 
[e.g., Buckmaster and Matalon (1988); Dold (1989, 1991)]. 

The initial presence of a vortex alters the picture considerably [Macaraeg, Jackson and Hussaini 
(1992)]. First, a hot spot develops within the viscous core of the vortex and evolves into an almost circular 
flame which grows with time according to a similarity rule and so does the reactant-consumption rate, as in 
the case of Marble. Next, the Linan and Crespo scenario or the tribrachial flame configuration must evolve 
at infinite distance to the right and left of the origin. The diffusion flames emanating from this configuration 
then moves toward and finally merges with the flame in the reacted core region. The above picture is only 
valid for the case when the adiabatic flame temperature is greater than either of the two initial temperatures. 
When the adiabatic flame temperature lies between the initial temperatures, a well defined ignition point 
does not occur, but rather the premixed region merges smoothly into the diffusion flame region, independent 
of the vortex Reynolds number. This scenario is consistent with the Linan and Crespo analysis. 

The presence of an acoustic field or time-dependent pressure perturbation affects both the ignition and 
the flame structure. To investigate the ignition regime, an asymptotic theory is developed here which holds 
strictly for near-equal temperatures of the reactants. For low frequencies of a time periodic pressure pertur- 
bation, ignition can be accelerated or significantly delayed depending on the phase. For moderate to high 
frequencies ignition is always enhanced. An approximation based on homogeneous ignition theory is given 
which predict these trends very well, and therefore is proposed as a good engineering approximation for 
obtaining the ignition times for the more complicated flow fields. For the complete time evolution from 
nearly frozen flow to near equilibrium flow, an asymptotic analysis based on small heat release is employed. 
In this limit, direct numerical simulations show that acoustics have little effect on the flow field in the vor- 
tex core, while the flow outside the core oscillates. It is believed that viscosity in the reacted core damps out 
the acoustic field. 

In the next section the problem is formulated. Section 3 presents an analysis of the ignition regime, 
valid only for near equal initial temperatures of the reactants, using a combination of large activation energy 
(and, hence, large Zeldovich number) asymptotics and numerics for several different pressure profiles. The 
pressure profiles chosen are either a linear, pulse, or a sinusoidal function. This last case is stressed since it 
represents a single Fourier component of a more general pressure disturbance. Also in this section, the 
approximation based on homogeneous ignition theory is presented. Section 4 presents selected numerical 
results of the continuous evolution from nearly frozen flow to near equilibrium flow for the sinusoidal pres- 
sure profile under the assumption of small heat release. A description of the diffusion flame regime is given 
in Section 5 and our conclusions are presented in Section 6. 

2 . PROBLEM FORMULATION 

In this section the problem for the time evolution of initially unmixed species occupying adjacent 
half-planes which are then allowed to mix and react in the presence of a vortex subject to pressure 
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disturbances is formulated. The nondimensional equations governing this field, assuming constant viscous 
and thermodynamic properties, are given by Buckmaster and Ludford (1982) and Williams (1985), 


p T = 1 + 5P c (t) + yM 7 P(x,y,z,l), (2.1a) 

P, + (pH), + (pv), + (pw) 2 = 0, (2.1b) 

p [ u, + u u x + v u y + w u 2 ] + P x = Sc ( V 2 « + y K x ), (2.1c) 

p [v, + u v x + v v y + w v z ] + Py = Sc (V 2 v + y K y ), (2. Id) 

p[w, + u w x + v Wy + w w 2 ] + Pi = Sc (V 2 w + y K z ), (2.1e) 


y — 1 dP c 

p[T, + uT x +vT } +wrj-5 JL ^ — 

- (Y- l)M 2 [P, + uP x + vPy + wP t ] = v 2 r + (y- l)M 2 5c d> + ptl, (2.10 

p [Fjj + u Fj' X + v F j<y + w Fj' Z ] = V 2 Fj - £2, j = 1,2, (2.1g) 

£2 = Da pF l F 2 e~ ZelT , (2.1h) 

where K = u x + v y + w 2 , <b is the viscous dissipation term, and V 2 is the three-dimensional Laplacian 
operator. Here, («,v,w) are the velocity components in the (x,y,z) directions, respectively; p is the den- 
sity; T is the temperature; and F { and F 2 the mass fractions of the fuel and oxidizer, respectively. The total 
pressure P = 1 + bP c + yM 2 P is written as the sum of the compressible component P c plus incompressi- 
ble component P, as is consistent with the small Mach number approximation [e.g., Majda (1984)], with 5 
the amplitude of the compressible component. The actual size of 5 will be chosen in the course of the 
analysis. The chemical model is assumed to be a one-step, irreversible Arrhenius reaction. The nondimen- 
sional parameters appearing above are the Schmidt number Sc = v / D assumed equal for both species with 
D the species diffusion coefficient and v the kinematic viscosity; the Zeldovich number Ze = E / (R° T^) 
with E the dimensional activation energy and R° the universal gas constant; the Damkohler number Da 
defined as the ratio of the characteristic time scale t d to the characteristic reaction time scale t R ; y the ratio 
of specific heats; and finally P the heat release per unit mass of F i >oa . In specifying the Damkohler number 
in the course of the analysis, the appropriate choice for the characteristic time scale is made. The density, 
temperature, pressure and mass fractions were nondimensionalized by their initial values 7^, 
P„ = P' M R° T„ and F l oo , respectively, that would prevail if there were no pressure disturbances. The 
Lewis number was assumed to be unity, which implies that Sc = Pr , where Pr = p„ C p v / X is the Prandtl 
number, X is the thermal conductivity and C p is the specific heat at constant pressure. Lengths and veloci- 
ties are referred to the relevant diffusion characteristic scales Q = D t d and U d ~l d ! t dy respectively. In 
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this case, the Reynolds number based on this choice of l d and U d is given by Re = l d U d /v = 1/Sc. The 
Mach number M is defined as U d ia 00 , the ratio of the characteristic diffusion speed to the speed of sound. 
The above choices for t d , l dy and U d are consistent with the scales chosen in the absence of a vortex. 


An analysis of the ignition zone with acoustic interactions must take into account the two fundamental 

ratios 

characteristic time , t d l d I U d 

x = = 

acoustic time , t a l a ld^ 


characteristic acoustic length 
diffusion length 


la_ 
Id ’ 


These two ratios are related via 


1 

N M ’ 


where M is the Mach number defined above. Here, our analysis is based on the assumption that t = 0(1), 
and hence the acoustic wavelength is much longer than the characteristic diffusion length scale for small 
Mach numbers. Thus, the pressure gradient within the combustion zone is negligibly small and the pressure 
is essentially a function of time, and so the flame is treated as isobaric to an excellent approximation. In 
addition, our analysis further assumes that the magnitude of the acoustic time relative to the reaction time is 
given by 

characteristic time , t d t a 

Da = = x — . 

reaction time , t R t R 


Taking x = 0(1), we shall assume that the reaction time is much faster than the acoustic time, and so we 
take the Damkohler number to be proportional to 


where Ze » 1 is the Zeldovich number defined above, thus defining the characteristic time. The exact form 
is chosen in the course of the analysis. In our analysis, if the Mach number M was not chosen small, we 
would also need to consider spatial variations of the pressure disturbance inside the reaction zone. 

McIntosh (1986, 1989, 1991a,b, and the references cited therein) pointed out the importance of the 
two fundamental ratios given above when the characteristic velocity is given by the flame speed Uf. For 
this case, three distinct cases of flame-acoustic interactions arise depending on the magnitude of x, and these 
have been classified in McIntosh (1991a). Further phenomena can exist depending on the magnitude of the 
acoustic time/reaction time ratio, and a discussion of this can be found in Clarke (1985). For x = 0(1), 
McIntosh and Wilce (1991b) point out that audible sound and its interference with flames would come under 
this category, and a considerable body of work exists using this assumption (e.g., references cited therein). 

We shall now assume that the flow is independent of the downstream direction x for all time. 
Although not true for any real flow situation, it is consistent with experimental observations that streamwise 


- 5 - 


vortical structures in turbulent mixing layers are essentially aligned with the flow and are of great extent. 
Also, neglecting the x -dependency can be thought of as a local approximation for the vortices found in tur- 
bulent flames. This type of approximation has also been used previously in describing the evolution of a 
shear flow with an embedded streamwise vortex [Corcos (1988), Pearson and Abernathy (1984)]. We shall 
also assume that the Mach number is small, as is typical for flames [e.g., Buckmaster and Ludford (1982), 
Williams (1985)]. Thus, the governing equations (2.1), in the limit of small Mach number and neglecting 

all dependency in x , reduce to 


pr = l + 5MO, 


(2.2a) 

P, + (pV)y + (PW^ =0, 


(2.2b) 

p [ u, + v u y + w u z ] = Sc V 2 h 


(2.2c) 

p [ v, + v v y + w v z ] + Py = 5c ( VS- + — K y ), 


(2.2d) 

p [ w, +vw y + w w 2 ] + P, = Sc ( V 2 w + y K z ), 


(2.2e) 

p[T, +vT y +wT z ]-d 1 ^- d( =v 2 r + pa, 


(2.20 

p [F Jtt + V Fj,y + w Fj^ ] = V 2 F; - £2, j = 1,2, 


(2.2g) 

where V 2 = ( ) yy + ( )„ is the two-dimensional Laplacian operator, K =v y + w z , and we have considered 
the distinguished limit M 2 « 8 « 1. Since we are interested in reacting flows with vortical motion, we 
transform (2.2) into cylindrical coordinates (y,z) -> (r,0) by using the transformations y = rcos0 and 
z = r sin 0. If we define g to be the radial velocity and A to be the tangential velocity, i.e., 
v = g cos 0 - h sin 0 and w = g sin 0 + h cos 0, then (2.2) transforms to 

pr = l + 8P e (f). 


(2.3a) 

P( + (P«)r + y(p /l )e+ -^ = °’ 


(2.3b) 

h 

p[u, + g u r + — U 6 ] = Sc V 2 !/ , 


(2.3c) 
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Plk + g h r + ~^e + ] + — Pq = 


r r 


Sc 




r 2 r 2 


3 r 


(2.3e) 


' IT, + g T r + A r j _ g 1_L <1 = V 2 7 


y dt 


+ pn. 


(2.30 


P ( F y.< + £ F j,r + ~ F j,e] = V 2 F 7 - - n, j = 1,2, (2.3g) 

where V 2 is now the two-dimensional Laplacian operator in cylindrical coordinates, and 

K = g r + g /r + h e /r. The system (2.3) must be solved subject to appropriate boundary and initial condi- 
tions. 


Since it is known that small pressure disturbances of 0{Ze~ l ) can cause 0(1) changes in certain flame 
properties such as the burning rate, extinction, transition from laminar to turbulence, etc. [Buckmaster and 
Ludford (1982), Buckmaster (1992), Kapila (1992)], we shall limit our attention to small pressure distur- 
bances by defining for the total pressure 

P = l + Ze ~ l ^r L J F c(t) + 0(M 2 ), 8 = — (2 4) 

the factor y l(y- 1) in the definition of 5 was chosen for convenience. Thus, we are concerned with 
O ( Ze ) acoustic perturbations superimposed on a basic state. The pressure perturbation P c (t ) can be deter- 
mined by solving the acoustic equations in the outer field [e.g., Van Harten, Kapila and Matkowsky (1984)], 
and indicates that the pressure can be controlled, for example, from the surroundings. A general expression 
for P c can be found in Majda (1984) and will depend on the type of conditions imposed on the geometry 
under consideration. Since we are interested in the influence of pressure variations on the flame, we will 
not solve this outer problem and, henceforth, consider the pressure disturbance P c to be prescribed. 

We note here that it is the density which couples the momentum equations to the temperature and 
mass fraction equations. As such, the numerical solution to the full system (2.3) is a formidable task. One 
approximation which will simplify (2.3) to a more tractable system is the assumption of small heat release 
[e.g., Van Harten, Kapila and Matkowsky (1984)]. The resulting system can then be solved numerically, and 
will be presented in Section 4. For the ignition problem, the assumption of small heat release is not needed 
in order to simplify the system, provided we restrict the analysis to near equal initial temperatures, and 
hence, densities, of the two gases. This analysis is presented in the following section. 

3. IGNITION 

At time t - 0, the reaction rate is exactly zero owing to the product F { F 2 = 0. For t > 0, the fuel and 
oxidizer begin to mix by diffusion, as well as by convection due to the presence of the vortex, and the reac- 
tion rate is no longer zero. For small time, it can thus reasonably be assumed that the effect of the reaction 
on the overall flow field is small. The solution for small time in which the reaction rate term is neglected is 
known as the inert or chemically frozen solution, and will be denoted by the superscript I. In what follows, 
we only consider the case of near equal initial temperatures. This assumption leads to the inert solutions 
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T 1 = p' = 1 + 0(Ze~ l ), where the 0(Ze~ l ) term is included to allow for small initial temperature 
differences that might exist in the absence of any pressure disturbances. Thus, the leading-order inert solu- 
tion corresponds to the constant density approximation prevalent in the combustion literature. As time 
increases, more of the combustible mixes until, at some finite time, a thermal explosion occurs characterized 
by significant departure from the inert. To analyze the ignition process, we determine the effect of the grow- 


ing reaction rate by expanding about the inert solution as 
7 = 1 + Ze- l [T x + P c (t)] + 0 (Ze~ 2 ), p = 1 + Ze~ l p, + 0(Ze~ 2 ), Fj = F j + 0(Ze~ l ), (3.1) 

u = u‘ + 0(Ze~ [ ), g=g' +0(Ze~ 1 ), h = h' + 0(Ze~ l ), (3.2) 

and take the asymptotic limit Ze -> The leading-order equations are given by 

Pi + T^^-jPAO, @.3a) 

gi+-h‘ B + ^-= 0, (3.3b) 

r r 


«/ + g 1 «/ + 7 - “e = Sc VV , 


(3.3c) 


g! + g' gr + ~ g9 ~ -y- + P r = ^ 




r 2 r 2 


(3.3d) 


h! + g 1 hi +—h' Q + + -7 e = Sc 


r r 


v ! /,'-4+4s’« 

r l r l 


(3.3e) 


h' 


T Ut + g‘ Tur + 7 - T i,e = V 2 ^ + A F\ F' 2 e 


/ rl Ji +p c«) 


(3.30 


J,l 


+ g ! F 1 


h< 


j- r 


7,0 


V 2 f', j = 1,2, 


where we have chosen the Damkohler number Da to be 


Da 


A 

pze 


.Ze 


(3.3g) 


(3.4) 


and A is some constant that will be chosen in the course of the analysis. This particular choice of the Dam- 
kohler number ensures that a distinguished limit exists, in that the reaction rate term is of the same order in 
a Zeldovich number expansion as the time derivative terms [see, e.g., Buckmaster and Ludford (1982)]. 
Note that the pressure disturbance now appears in the exponential and thus, will play a crucial role in the 
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ignition process. 

We first see that the continuity and momentum equations (3.3b-e) are decoupled from the rest of the 
system and thus, will completely determine u ! , g 1 , h 1 and P . At ( = 0, we situate a point vortex at the ori- 
gin and allow it to diffuse under the action of viscosity. Then, assuming the flow in the (r,0) plane to be 
axisymmetric, the induced flow field is given by [Lamb, 1932] 

g' = 0, h* = [1 — cxp( — r^/ 4 1 Sc )] , (3.5) 

P(r,t) = - J — ( ' ,f) dr, (3.6) 

r r 

and u 1 is found from the equation 


«/ + — «e = Sc V 2 w ; . 


(3.7) 


The solution (3.5)-(3.6) is the solution for the incompressible Oseen vortex, where R = T/2kv is the vortex 
Reynolds number and T is the circulation parameter. The vortex Reynolds number R is considered as a free 
parameter and its choice determines the importance of the vortex in the initial mixing of the reactants. 
When R = 0 (i.e., T = 0), there is no vortex and the initially separated reactants mix by diffusion only. 
When the vortex Reynolds number R is non-zero, the initially separated reactants mix by diffusion and con- 
vection. The solution to (3.7) has been discussed previously in the context of streamwise-spanwise vorticity 
interactions for nonreacting flows [Corcos (1988), Pearson and Abernathy (1984)], Note that u 1 does not 
appear in the temperature or mass fraction equations since we neglected variations in x coupled with the 
zero Mach number approximation, and so the analysis is valid for either streamwise or spanwise vortices. 
The velocity component u 1 will influence the temperature equation if true compressibility effects are taken 
into account; such will necessarily occur after ignition has taken place. 

With g 1 and h ] now known, the following equations can be solved for T { and Fj 




T\ a = V 2 T\ + A F[ F' 2 e T ' +PcU \ 


(3.8) 


Fj,i + ~Fj,e = V 2 Fj, 7 = 1,2. (3.9) 

The appropriate boundary and initial conditions are given by 
T l = 0, F\ = 1, F *2 = 0 at t = 0, r > 0, 0 < 0 < 7t, and t > 0, r oo, 0 < 0 < 7t, 

T\ ~ Pr> F\ = 0, F 2 = s <(> at t = 0, r > 0, n < 0 < 2n , 


and t > 0, r — > oo, % < 0 < 2 n. 


(3.10) 
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where § is the equivalence ratio defined as the ratio of the initial mass fraction of the fuel F u „ to the initial 
mass fraction of the oxidizer F 2 ,-»> and is the parameter which allows for small initial temperature 
differences in the absence of pressure disturbances. If <)> = 1» th e mixture is said to be stoichiometric; if 
<1> > 1 it is fuel rich; and if <|> < 1, it is fuel lean. Also, if (V is less than zero, the oxidizer (species 2) is 
relatively cold compared to the fuel (species 1); and if [3/ is greater than zero, it is relatively hot. As t 
increases, the solution for T { becomes unbounded at some finite time (t , s ) and location (y ig , z ig ). This 
characterizes the ignition regime. The special case R = 0, in which two initially unmixed species are 
allowed to diffuse without the mixing generated by the vortex, will be discussed in subsection 3.1, while the 
case for R >0 will be presented in 3.2. Finally, with T x determined, the density perturbation pj can be 
found from (3.3a). 

3.1. R = 0. In the absence of any vortical motion, the governing equations reduce to the 
reactive-diffusion equations 


t -t , a f' f 1 + 

* l,t — * \,zz + ™ r i r 2 € > 




= F 


J ,zz * 


j = 1,2, 


subject to the boundary and initial conditions 

r, = o, F\ = 1, F' 2 = 0 at t = 0, z > 0, and t > 0, z 

r, = p T , F[ - 0, F' 2 = (J> _1 at t = 0, z < 0, and l > 0, i -» 

The parabolic equations for the mass fractions can be integrated to yield 

F\ = 'F, F' 2 = r 1 (1 - TO, 


(3.11) 

(3.12) 


(3.13a) 

(3.13b) 


(3.14) 


where 


T 7 = — (! + er/ q ), 


(3.15) 


and q is the similarity variable, defined as 


q = 


2<I • 


(3.16) 


Upon substituting the inert solutions into the temperature perturbation equation, and transforming to the 
similarity plane, a single equation for T i emerges, given by 


4/r M -r Mt ,-2qr^ = 4t'F(l-'F)e 


r,+P c (0 


(3.17) 


which must be solved numerically subject to the boundary and initial conditions 
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Ti= 0 at t = 0, z > 0, and T > 0, z -» <*>, (3.18a) 

Tj = p 7 a/ t = 0, z < 0, and T > 0, z -°o (3.18b) 

Here, the constant A appearing in the definition of the Damkohler number has been chosen to be <|> for con- 
venience. The special case of no pressure disturbance (P c = 0) corresponds to the results of Linan and 
Crespo (1976) which has been discussed in the Introduction. Thus, the results of this subsection can be 
thought of as an extension of the Linan and Crespo problem to include pressure disturbances. In the follow- 
ing subsections we will select three different pressure variations, and compare with the results of Linan and 
Crespo. The first two choices for the pressure variation presented below corresponds to the choices taken by 
Ledder and Kapila (1991) in their study of the response of premixed flames to pressure perturbations. The 
last choice will be a sinusoidal pressure variation, which represents a single Fourier component of a more 
general pressure disturbance. 

3.1.1. LINEAR PRESSURE DECREASE Let the flow be subject to a linear pressure decrease at 
t = 0 of the form 


/% (/ ) = - co r , (3.19) 

where 1/co > 0 charactererizes the time scale of the pressure disturbance. The influence of this linear pres- 
sure profile on premixed flames has been discussed by Ledder and Kapila (1991), in which extinction was 
seen to occur for all values of co. This profile was chosen for this study because of its simplicity in form in 
that the ignition time and location will only be a function of the parameters co and p r , and also may 
represent the first term in a Taylor Series expansion of a more general pressure disturbance. Figure 1 is a 
plot of the maximum of T j versus / for several values of co and p r = 0. The case co = 0 corresponds to that 
of Linan and Crespo (1976) in which no pressure disturbance is present. As co is increased from zero, the 
ignition time also increases until, at some finite value co s co c = 0.1395, the ignition time is pushed to 
infinity. Thus, small values of co can significantly delay ignition. For values of co greater than to c , the 
maximum of T { initially increases, reaches a maximum, then begins to decrease back to zero as time 
increases further. Thus, ignition does not occur in the classical sense. Mathematically, when the pressure 
disturbance becomes large and negative, the nonlinear source term of (3.17) becomes exponentially small 
and so can be neglected; T x then evolves according to a homogeneous heat equation. Although the solution 
of this system suggests that ignition will not occur for co greater than co c , we note that the presence of an 
infinite amount of fuel and oxidizer necessitate that ignition will eventually take place. Asymptotically, the 
expansion (2.4) breaks down when t is of the order 0(Ze ), and a new expansion must be employed, one 
which allows 0( 1) changes in the pressure and the temperature; the waves now cease to be acoustic in 
nature. The leading-order system in this case is no longer that with constant density, and hence, the full sys- 
tem (2.3) must be solved numerically. This numerical problem is beyond the scope of the present study. 
Figure 2 is a plot of the ignition time t lg versus co for three values of the temperature ratio parameter p r . 
The case of p r = 0 was presented in a slightly different manner in Figure 1. Here, it is clear that as co 
approaches 0.1395 from below, the ignition time goes to infinity, and for co greater than 0.1395 there is no 
ignition. The case to = 0 corresponds to that of Linan and Crespo (1974); ignition time increases as (3 r is 
decreased from zero (F 2 ,« is relatively cold compared to F x and decreases as p r is increased from zero. 
For each value of the temperature ratio parameter, a critical value of co is seen to exist; this value being less 
than 0.1395 for negative p r and greater than 0.1395 for positive values of p T . In the absence of any 
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pressure disturbance, Linan and Crespo have shown that ignition takes place along a line parallel to the y- 
axis which is located at z ig = 0 for (3 r = 0 and resides in the hotter region for ^ 0. In the presence of the 

linear pressure disturbance and for p r = 0, ignition always occurs along the line z ig = 0 for any value of co. 
However, for p T * 0, the ignition line z = z lg (z ig < 0 for p T > 0, and z ig > 0 for p r < 0) moves toward the 
line z = 0 as co is increased. 

3.1.2. PRESSURE PULSE Let the flow be subject to a pressure pulse at t = 0 of the form 

P c (t) = 4 P A (2"“' 2 - 2 _2 “' 2 ), (3.20) 

where oT 1/2 > 0 is the pulse width and | P A | is the amplitude of the pressure pulse. This corresponds to a 
smooth increase beginning at zero, increasing or decreasing to an extremum | P A | at cor 2 = 1, and finally 
returning to zero. The influence of this pressure pulse on premixed flames has been discussed by Ledder 
and Kapila (1991), in which extinction is seen to occur for both negative and positive amplitudes, depending 
upon the parameters of the problem. Figure 3 is a plot of ignition time t lg versus co for various values of the 
amplitude P A with p r = 0. The case of P A = 0 corresponds to that of Linan and Crespo (1976) and is 
shown in the figure as a dashed line for reference. In addition, all data curves shown in this figure begin at 
t ig = 5.816 at co = 0, again consistent with the results of Linan and Crespo. For fixed P A < 0, the ignition 
time initially increases as co increases from zero, reaches a maximum, and then begins to decrease as co is 
increased further, reaching an asymptote as co °o. For fixed co, decreasing the pressure amplitude from 
zero increases the ignition time. Thus, ignition is delayed for negative pressures. In contrast to the previous 
case, for a positive pressure pulse, the ignition time decreases monotonically as co is increased from zero 
and asymptotes to a finite non-zero value as co — > <»; hence, ignition is enhanced. Since [J^ = 0, ignition 
takes place along a line parallel to the y-axis located at z ig = 0 for any value of P A and co. The general 
behaviour described in Figure 3 for p r = 0 is seen to hold for any value of p r , the only exception being the 
starting value of t lg at co = 0; for (3 r < 0, the starting value is above that given in the figure, while for 
Pr > 0 it is below. In addition, the location of the ignition line resides in the hotter region and moves 
toward z = 0 as co is increased from zero for fixed P A . 

An interesting feature is seen to occur as P A is further decreased below -2. This is shown in Figure 4 
where we plot t ig versus co for P A =-4 and p r =0. The ignition time increases slowly up to about 
co = 0.0018, where a sudden jump in t ig occurs. As co increases still further, the ignition time reaches a 
maximum and then begins to decay slowly, reaching an asymptote as co — » To better explain this jump 
phenomenon, we plot in Figure 5 the maximum of T l versus time for several values of co. Note that there 
is a sudden jump in the ignition time between co = 0.0018 and co - 0.002. For values of co below 0.0018, 
the maximum temperature perturbation is seen to increase monotonically, becoming infinite as ( approaches 
t ig . For values of co greater than about 0.0018, the maximum temperature perturbation increases at first, 
decreases slightly, and then increases sharply and becomes unbounded as t t ig . In the region where the 
maximum decreases, the pressure disturbance dominates in the exponent of (3.17), yielding a region where 
the chemical activity is suppressed, and so the temperature evolves according to a homogeneous heat equa- 
tion; this parallels that of the linear pressure case. However, unlike the linear pressure case, the pulse does 
not grow unbounded, but instead reaches a minimum and then increases back to zero. Thus, there exists 
some finite time at which the pressure no longer dominates the exponent, and the problem begins to once 
again evolve according to classical thermal explosion theory. A comparison of the pressure pulse and 
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maximum temperature perturbation is provided in Figure 6 for one selected value of co. Note that the max- 
imum in T j decreases over most of the pressure pulse, and that the temperature again rises once the pressure 
increases to some critical value. Finally, we note here that the above behaviour was seen to occur numeri- 
cally for all pressure amplitudes less than about -2.8, and also occurs for non-zero values of |3 r . 

3.1.3. SINUSOIDAL PRESSURE VARIATION Let the flow be subject to a sinusoidal pressure 
variation at t = 0 of the form 


P c ( 0 = P a sin (to t + O), (3.21) 

where to > 0 is the frequency, P A is the amplitude, and O is the phase of the pressure variation. This profile 
was chosen since it can be regarded as a single Fourier mode of a more general pressure disturbance. Figure 
7(a,b,c) are plots of t ig versus the phase angle d> for P A = 1, 2, 4, respectively, and for various values of to 
with p r = 0. The case of P A = 0 corresponds to that of Linan and Crespo (1976) and is shown in each 
figure as a dashed line for reference. We begin by noting that there are several important features which 
are common to all three figures. These are: (i) the maximum ignition time at a fixed amplitude corresponds 
to zero frequency with a phase angle greater than n (corresponding to an initial negative pressure wave), 
and that this maximum decreases and shifts to lower phase angles as co is increased from zero; (ii) for fre- 
quencies co - 1 and greater, the ignition time always lies below that of the dashed line for any phase angle 
and amplitude, indicating that ignition is enhanced over the zero pressure disturbance case; and (iii) for the 
high frequency range (co > 5), the ignition time is approximately a constant over the entire range of phase 
angles, and that this constant value decreases as the pressure amplitude P A increases. 

The same interesting feature that occurred for the pressure pulse case is also seen to occur for the 
sinusoidal pressure case, in that the ignition time undergoes a sudden jump in its value for certain combina- 
tions of amplitude, frequency, and phase. This jump phenomenon in the ignition time is clearly shown in 
Figure 7c for P A = 4, For the curves co > 0.1, the ignition time drops slowly as is increased from zero, 
reaches a minimum, and then begins to increase slowly as O is increased further until a critical value of 
is reached, where the ignition time then undergoes a sudden jump in its value. As 0 increases still further, 
the ignition time reaches a maximum and then begins to decrease back to its original value at O = 0. To 
better explain this jump phenomenon, we plot in Figure 8 the maximum of T { versus time for P A =4, 
co = 0.1 and O = it. The corresponding pressure wave is displayed as the dashed curve. In the region where 
the maximum temperature perturbation is almost constant, the pressure disturbance dominates in the 
exponent of (3.17) yielding a region where the chemical activity is suppressed, and so the temperature 
evolves according to a homogeneous heat equation, and ignition is delayed. As the pressure rises to some 
critical value, chemical activity resumes and ignition soon takes place according to classical thermal explo- 
sion theory. 

Figures 9 and 10 show that regions of chemical inactivity can exist in the ignition regime even in the 
high frequency limit. Figure 9 shows the maximum of T x as a function of time for P A = 2, co = 5, and 
<X> = it/ 10, while Figure 10 is a similar graph except now co = 10. In each graph, the corresponding pres- 
sure variation is shown as a dashed curve. Note that there are now several regions where the maximum 
temperature perturbation decreases and hence chemical activity is suppressed, with the pressure variation 
being negative in each of those regions. 
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Finally, we note here that the general behaviours described in Figures 7-10 for (3 r = 0 are seen to 
occur for any value of p 7 . The only exception being for (3 7 < 0, the ignition time is greater than that given 
for p = 0, while for p 7 > 0 it is less. 


3.2. R > 0. When vortical motion is taken into account, the governing reactive convection- 
diffusion equations for the temperature and mass fraction perturbations are 


Ti, ( + -exp(— r 2 /4t5c)]r lj0 = V 2 ?! + (t>F ; , f' 2 e 1 c °. 


(3.22) 


R Sc 


i.< 


[1 - exp(-r 2 /4f 5c)]Fy_ e = V 2 Fy, j - 1,2, 


(3.23) 


subject to the conditions (3.10). This system was solved numerically for a range of R , Pr = 0, Sc = 4> = 1. 
and with an assumed form for the pressure variation term P c {t). Implementation of boundary conditions is 
facilitated if the system is recast in Cartesian coordinates. The solution technique is a 2nd-order finite 
difference scheme on a nonuniform mesh. To resolve the structure in the core region of the field, a coordi- 
nate stretching is used. To avoid the singularity at the origin, no mesh points are placed there. The outer 
boundaries are set at 50 or 200 in the y -direction, and 20 or 50 in the z -direction. Grid resolution studies 
which at least doubled the computational mesh were carried out to ensure that structures were well resolved. 
The resolutions required ranged from a 64 2 mesh to a 256 2 mesh for large vortex Reynolds number. The 
time-stepping scheme is a four-stage Runge-Kutta which is formally 2nd-order but has an extended stability 
region making it accurate and robust for moderately stiff problems. All runs were performed on a Cray 
YMP. 

As t increases, the solution for T x becomes unbounded at some finite time (t lg ) and location (y ig , z lg ). 
This characterizes the ignition regime. The special case R = 0, in which two initially unmixed species are 
allowed to diffuse without the mixing generated by the vortex, has been presented in subsection 3.1. For this 
case, ignition takes place along a line parallel to the y -axis which is located at z, s = 0 for p 7 — 0 and 
resides in the hotter region for p r * 0. For R > 0, previous work [Macaraeg, Jackson, and Hussaini (1992)] 
has shown, in the absence of a pressure disturbance, that ignition occurs at a point rather than along an 
entire line. The ignition point is located at the origin for p 7 = 0 and for any value of R , and resides in the 
hotter region for p 7 * 0. For p 7 * 0, the ignition location spirals clockwise towards the viscous core center 
as R is increased from zero, with ignition taking place within the core for R > 70. 

To determine the combined effects of the vortex and a pressure disturbance on the ignition time, we 
produce here selected results using the sinusoidal pressure function (3.21) for P c (t ). In all calculations, we 
fix Sc = <j) = 1 and p 7 = 0, and vary the parameters P A , R , to, and d>. Since Pr = 0 and R > 0, ignition 
will occur in the vortex center. Figure ll(a,b) are plots of t ig versus the phase angle <b for P A - 2 and 4, 
respectively, and for vortex Reynolds number R — 0 (solid) and R = 100 (dash) with co = 0.5. The results 
with R = 0 (solid) are taken from Figure 7 and given here as a reference when R > 0. From Figure 1 1 we 
see that the effect of increasing R decreases slightly the ignition time for any phase. This is consistent with 
previous results obtained when pressure disturbances was ignored [Macaraeg, Jackson, and Hussaini (1992)]. 
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To illustrate the structure in the high frequency limit, we plot in Figures 12 and 13 the maximum of 
the temperature perturbation T x as a function of time for to = 5 and 10, respectively, and for various values 
of the vortex Reynolds number R with P A = 2, = it! 10, and p r = 0. The special case of R = 0 was pre- 

viously shown in Figures 9 and 10, respectively, and is used here as a reference case for when R > 0. In 
both figures, the structure is similar to the R =0 case, except that increasing R decreases the ignition time 
slightly until about R = 100, when the ignition time remains constant as R is increased further. As stated 
earlier, the temperature maximum always occurs in the vortex center. 


3.3. APPROXIMATION FOR THE IGNITION TIME. Since a certain amount of mixing must 
take place before ignition, a surprisingly good approximation to the ignition time in the presence of a pres- 
sure disturbance can be found by considering the "homogeneous" problem 


dT 1 1 T, +/>.(/) 

— p i c 

dt t ° 


(3.24) 


where the factor is the ignition time in the absence of a pressure disturbance chosen to be the Linan and 
Crespo (1976) result, which accounts for the role that diffusion plays in the mixing of the reactants as well 
as the role that the initial temperature profile plays. For p r = -2, 0, and 2, we find that = 13.1961, 
5.8285, and 1.7736, respectively, for our choice of the time scales (3.4). Integrating (3.24), we find that 


T l = -In 




p c 0) 
e dt 


(3.25) 


l ig 0 


Ignition takes place when the quantity in the brackets of (3.25) vanishes; that is, the "homogeneous" ignition 
time is given implicitly by the equation 


j/ e<W dT = /£. (3.26) 

0 

In the absence of a pressure disturbance, P c (t) = 0 and therefore t" = t°. Below we give selected com- 
parisons between the ignition times predicted by this "homogeneous" theory and by the full numerical solu- 
tions of the previous subsections. 

First, consider the simple linear pressure profile given by (3.19). Substitution into (3.26) yields for the 
approximate ignition time 




(3.27) 


Note that (3.27), in addition to defining the approximate ignition time, also approximates the critical value 
of to for which ignition can take place. That is, ignition occurs for 0 < to < 0)^, where co c w = 1 / tR is the 
critical value; for co > <aR, ignition does not take place. In Figure 14(a) we plot the ignition times for the 
numerical solutions of subsection 3.1.1 (solid) and the approximate solutions (dashed) for (3 r = -2, 0, and 2 
as a function of co. For each p r , note that the approximate ignition times compare rather well to that of the 


- 15 - 


full solutions. Now consider the sinusoidal pressure function given by (3.21). For this case, (3.26) must be 
solved implicitly for /£ for fixed values of P A , CO, <& and p r . In Figures 14(b,c) we plot the ignition times 
for the numerical solutions of subsection 3.1.3 (solid) and the approximate solutions (dashed) as a function 
of phase and for P A = 4, to = 0.1 and P A = 4, co = 0.5, respectively, with p T = 0. In each figure note how 
well the approximate ignition times compare to the ignition times of the full numerical problem. Thus, we 
see that our theory of "homogeneous" ignition predicts surprisingly well the true ignition times. Finally, we 
note here that our definition of the "homogeneous" ignition time can also be used to predict ignition when 
vortical motion is present, since ignition times vary little with vortex Reynolds number R . 

4. DIRECT NUMERICAL SIMULATION 

To investigate post-ignition events up to the development of a flame sheet, the direct numerical solu- 
tion to (2.3) is necessary. Unfortunately, as mentioned above, the numerical solution to the full system (2.3) 
presents a formidable task. One approximation which will simplify (2.3) to a more tractable system and has 
often appeared in the literature is the assumption of small heat release [e.g., Matkowsky and Sivashinsky 
(1979); Van Harten, Kapila and Matkowsky (1984)]. The resulting system can then be solved numerically. 
We begin the analysis of small heat release by assuming P « 1. The temperature and density must now be 
rescaled according to 

r-> 1 + P(r + P c (t)), p-*l + pp, ( 4 - ] ) 

while the other variables remain unchanged. To leading-order in p, the system (2.3) reduces to 

T, + -^| £ [ l-exp(-r 2 /4tSc)]r 0 = V 2 T + <l)r 1 F 1 F 2 e X(7 ’ + Pc( ' )) , (4.2) 

r 2 


Fj.t + 


R Sc 


[1 - exp(-r 2 /4/ Sc)]Fj t Q = V 2 F ; - §X 1 F [F 2 e 


xcr + /yr» 


j = 1,2, 


(4.3) 


for the temperature and mass fractions, while the leading order solution to the continuity and momentum 
equations (2.3b)-(2.3e) is the incompressible vortex solution (3.5)-(3.7). The density can be found from the 
gas law. In the above equations, X= pZ<? is the product of the heat release parameter and Zeldovich 
number, 5 = py/(Y“ 1), and (3.4) with A = <|> has been used in the definition of the Damkohler number. 
Finally, we note here that by using the definition (4.1) for the temperature, we have subtracted out the time 
varying part due to pressure fluctuations; that is, if T oscillates at any point in the flow field, it will be due 
solely to the nonlinear source term in (4.2). 

The system (4.2)-(4.3) must be solved numerically subject to the boundary and initial conditions 

7 = 0, F 1 = 1,F 2 = 0 at t = 0, r > 0, 0 < 0 < 7T, and t > 0, r -> <», 0 < 0 < n, (4.4a) 

T = p r , F x = 0, F 2 = <t> -1 at ^ = 0, r > 0, < 0 < 2 tc , and t > 0, r -> «>, n < 0 < 2tt. (4.4b) 

This system was solved by the finite difference scheme mentioned in Section 3. To illustrate the numerical 
solution of the continuous evolution from nearly frozen flow to near equilibrium flow, we produce here 
selected results using the sinusoidal pressure function (3.21) for P c (t). In all calculations we take 
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Sc = <|> = 1 and p r = 0. 


To investigate the combined effects of a pressure disturbance and of vortical motion on the structure 
of the full solution, we begin by examining the structure in the absence of these effects. Figure 15(a) is a 
plot of the maximum values of T as a function of time for R = 0, P c (t) = 0 and X = 1, 0.1 and 0.01. Since 
Pr = 0 and R =0, this maximum occurs along the line z = 0. For X = 1, the temperature increases gradu- 
ally from zero at t = 0 to T - 0.5 as t — > this maximum value corresponding to the flame sheet value 
given in the next section. As X decreases from one, the rise in temperature near t = 0 increases sharply, 
indicating the existence of a well-defined ignition regime. An interesting trend develops if pressure fluctua- 
tions are now included. In Figure 15(b) we plot the maximum values of the temperature as a function of 
time for R = 0 and for the sinusoidal pressure function (3.21) with P A =4, co = .5 and <I> = 1 .25 tc. For 
X = 1 and 0.1, the temperature rises from t = 0 to a maximum value slightly below the flame sheet value of 
0.5, and then begins to oscillate in time with decreasing amplitude; as / -> °°, the oscillations die out and 
the temperature maximum approaches the value 0.5. As X is decreased from 0.1 to 0.01, ignition takes 
place so rapidly that the acoustic field is seen to have almost no effect on the temperature field. The 
absence of any effect due to pressure disturbances is consistent with the flame sheet solution presented in 
the following section. A surprising phenomenon occurs when vortical motion is included. In Figure 15(c) 
we again plot the maximum values of the temperature as a function of time, again for X = 1, 0.1 and 0.01, 
with R = 50 and the same pressure profile as in Figure 15(b). Since p r = 0 and R > 0, this maximum 
occurs in the vortex center. Note that for all three values of X, the temperature profile evolves in time with 
little influence from the pressure field. That is, the temperature within the viscous core, once it reaches a 
maximum, remains almost constant with time, indicating that viscous effects damp out the acoustics. In con- 
trast, we show in Figure 16 the maximum temperature at the origin (solid) and the maximum temperature 
far away from the core (dashed) as a function of time for the same conditions as in Figure 15(c) with X = 1. 
From this figure we see that, although the temperature within the core remains almost constant with time, 
the outer flow is oscillating. This is consistent with Figure 15(b), since the outer flow does not feel the 
influence of the vortex so that the acoustic field has a noticeable effect. Finally, the oscillations in the outer 
flow will eventually die out, consistent with the flame sheet solution given in the next section. 

5. DIFFUSION FLAME REGIME 

After ignition has taken place, a thin diffusion flame exists and is characterized by a chemical reaction 
time that is much smaller than a characteristic diffusion time. Chemical reactions then occur in a narrow 
zone between the fuel and the oxidizer, where the concentrations of both reactants are very small. 
Mathematically, the assumption of very fast chemical reaction rates leads to the limit of infinite Damkohler 
number which reduces the diffusion flame to a flame sheet (i.e., local chemical equilibrium). This assump- 
tion significantly reduces the complexity of the problem since it eliminates the analysis associated with the 
chemical kinetics. For many flows, the assumption of local chemical equilibrium adequately predicts the 
location and the shape of the diffusion flame [Buckmaster and Ludford (1982); Williams (1985)]. For finite 
values of the Damkohler number, equations (4.2) must be solved numerically, and this was done in the pre- 
vious section. 

We begin the description of the diffusion flame regime by defining the following conserved variables 
[see, e.g., Williams (1985)]: 
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r+F, = p T +(l-p T )Z, 

(5.1) 

T + F 2 = (Pj- + <() *) (1 — Z), 

(5.2) 

„ f,-f 2 + 4>- 1 

^ — ■ r 

1 /K-l 

(5.3) 


where Z is the mixture mass fraction and satisfies the convection-diffusion equation 

Z, + -^[1 -exp(-r 2 /4f Sc)]Z e = V 2 Z, (5.4) 


subject to the initial and boundary conditions 

Z = 1 at t = 0, r > 0, 0 < 0 < n, and t > 0, r -> °°, 0 < 0 < n, (5.5a) 

Z=0 at t = 0, r > 0, n < 0 < 2rc, and t > 0, r °° n < 9 < 2jr. (5.5b) 

In the limit X -» 0 (i.e., infinite Damkohler number) the flame sheet solution is given by 

F , = 1 - (1 + ((T 1 ) (1 - Z), F 2 = 0, (5.6a) 

7 = (Pr +<t>-')(l-Z), (5.6b) 

valid for Z > Zf, and 

F , = 0, F 2 = <1)- 1 -(1-m))- 1 )Z, (5.7a) 

T = p r + (1 - Pr)Z, (5.7b) 


valid for Z < Zf, Here, Zf defines the location of the flame sheet where both the reactants vanish, given by 
the implicit relation 


Zf = 


1 

1 + <f>’ 


and T takes the adiabatic flame value 


Tf = 


l + Pj’ <j) 

! + <()' 


(5.8a) 


(5.8b) 


Note that the flame location is independent of P; and p. Once Z = Z(r,0, t) is known, then the other vari- 
ables ( T,F i,F 2 ) can be found from (5.6)-(5.7). We remark here that if p T = 0 and <p = 1, then T ( = 0.5, 
which is consistent with Figures 15-16. We also note that the flame sheet structure and position is 
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independent of the pressure disturbance P c (/), again consistent with Figures 15-16 in the limit i — > 

As mentioned in the Introduction, equation (5.4) for Z is exactly Marble’s problem [Marble (1985)] 
when ignition events are ignored (i.e., the vortex is turned on at t = 0 with the initial configuration given by 
the flame sheet), though derived for order one heat release and with the constant density approximation, and 
was solved numerically by Laverdant and Candel (1988) and Rehm, Baum, Lozier and Aronson (1989). 

6* CONCLUSIONS 

The problem of a flame interacting with a vortex is supposed to model certain fundamental mechan- 
isms of turbulent diffusion flames which form the basis of many propulsion devices. Acoustic waves are 
always present in these combustion systems, so that it is important to study the role of acoustics in such 
flame/vortex interactions. Having said this, we refrain from any further comment on the implications of the 
present study to practical applications. 

The present study assumes a one-step, irreversible Arrhenius reaction between initially unmixed 
species occupying adjacent half-planes which are allowed to mix and react by convection and diffusion in 
the presence or absence of a vortex, and investigates the influence of the acoustic field (spatially uniform 
but time dependent pressure waves of small amplitude) on the ignition time and flame structure which are 
functions of vortex Reynolds number, initial temperature difference of the reactants, equivalence ratio and 
Schmidt numbers. 

In the case of zero vortex Reynolds number and near equal initial temperatures of the reactants, this 
study concludes that the low frequency pressure waves accelerate or decelerate ignition depending on the 
phase; the pressure waves with moderate to high frequencies always enhance ignition and are essentially 
independent of phase. Another key result is that the simple theory of homogeneous ignition developed in 
this study is found to predict these trends very well, and therefore it is proposed as a good engineering 
approximation for obtaining the ignition times for the more complicated flow fields. When vortical motion 
was included, the overall ignition time decreases slightly and then asymptotes to a finite value with increas- 
ing vortex Reynolds number, consistent with previous results obtained by ignoring acoustics [Macaraeg, 
Jackson, and Hussaini (1992)]. While the asymptotic analysis is confined to incipient ignition, a direct 
numerical simulation is carried out to investigate post ignition events. This is limited by the small heat 
release assumption. Relaxing this rather drastically restrictive assumption is worth a serious attempt. An 
important conclusion from the numerical investigation is that acoustics has little effect on the flow field in 
the vortex core, while the flow outside the core oscillates. It is believed that viscosity in the reacted core 
damps out the acoustic field. 
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Figure 3. Plot of the ignition time t ig as a function of to for several values of P A and for p r = 0. The dash 
curve corresponds to the results of Linan and Crespo (1976). Pulse Pressure Profile. 
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Figure 4. Plot of the ignition time t ig as a function of co for P A - -A and for p r = 0. Pulse Pressure Profile. 



Figure 5. Plot of the maximum of the temperature perturbation T x as a function of time t for several values 
of o and for P A = -4 and = 0. Pulse Pressure Profile. 
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Figure 6. Plot of the maximum of the temperature perturbation Tj as a function of time for to - 0.005, 
P A ~-4, and p r = 0. The dashed curve is the pressure profile. Pulse Pressure Profile. 




PHASE 


Figure 7. Plot of the ignition time t, g as a function of phase <I> for several values of the frequency to and for 
p r = 0; 7(a): P A - 1; 7(b): P A = 2; 7(c): P A = 4. The dashed curve corresponds to the results of Linan 
and Crespo (1976). Sinusoidal Pressure Profile. 
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Figure 9. Plot of the maximum of the temperature perturbation T x as a function of lime for P A = 2, co = 5, 
<I> = 71/ 10, Pr = 0. The dashed curve is the corresponding sinusoidal pressure profile. 



Figure 10. Plot of the maximum of the temperature perturbation T x as a function of time for P A = 2, 
co = 10, O = nf 10, = 0. The dashed curve is the corresponding sinusoidal pressure profile. 
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Figure 12. Plot of the maximum temperature perturbation Ti as a function of time for P A = 2, (0-5, 
<D = it / 10, p T = 0 and for (a) R = 10. (b) R = 50, and (c) R = 100. The dashed curve is the corresponding 
sinusoidal pressure profile. 
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Figure 13. Plot of the maximum temperature perturbation T x as a function of time for P A = 2, co = 10, 
<X> = n/ 10, p r = 0 and for (a) R = 10, (b) R = 50, and (c) R = 100. The dashed curve is the corresponding 
sinusoidal pressure profile. 
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Figure 14. Plot of the ignition times corresponding to the full numerical solutions (solid) and the homo- 
geneous ignition times (dashed) as a function of (a) co using the linear pressure profile; (b) phase d> using 
the sinusoidal pressure profile with P A = 4 and to = 0.1; and (c) phase O using the sinusoidal pressure 
profile with P A = 4 and to = 0.5. 
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